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We extend one dimensional asymmetric simple exclusion process (ASEP) to a surface and show 
that the effect of transverse diffusion is to induce a continuous phase transition from a constant 
density phase to a maximal current phase as the forward transition probability p is tuned. The 
signature of the nonequilibrium transition is in the finite size effects near it. The results are compared 
with similar couplings operative only at the boundary. It is argued that the nature of the phases 
can be interpreted in terms of the modifications of boundary layers. 
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I. INTRODUCTION 

History has shown us that the study of model sys- 
tems or toy models of real physical systems is the first 
step towards a deeper understanding of working of real 
physical systems [lj. In this spirit, asymmetric simple ex- 
clusion process (ASEP) is a prototypical model of non- 
equilibrium statistical mechanics that deals with systems 
with currents flowing through them. Such systems are, 
in general externally driven, for example, in living cells, 
motor proteins, traffic flows, driven diffusive systems, 
transport in condensed matter and mesoscopic systems 
etc[|1l. 

ASEP is comprised of particles moving in a particu- 
lar direction with the constraint of no two particles at 
the same site at the same time, called simple exclusion. 
A particle can hop if the next site is empty. Particles 
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FIG. 1: (a) Schematic diagram of 2-D asymmetric simple 
exclusion process (ASEP), with forward jump probability p 
(thick lines) and transverse (dashed lines) excursion proba- 
bility q, with p + 2q — 1. For q = 0, one gets decoupled 1-D 
ASEP. (b)The phase diagram of I-D and 2-D exclusion pro- 
cess. M = maximal current phase (p = 0.5), in the Q-phase 
the bulk density is a, and in the 7-phase it is 7. The dotted 
(thin solid) line is the first order (continuous) phase bound- 
ary. For the 2-D or the modified f-D, the M region widens 
with p inducing a transition to M in the shaded region. In the 
a-phase region p c — 2a and in the 7-phase region p c = 2/3. 
The points marked refer to Fig. [3] 
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are fed at one end, say i = 1 at a rate a and with- 
drawn at i = m[m — > 00) at a rate (3 = 1 — 7 so that 
there is a current through the track. See Fig. [lja). 
The main interest in ASEP has been in the steady state 
properties, especially the nonequilibrium phase diagrams 
and the stability of phases as the external parameters 
or drives are changed. The phase diagrams in several 
cases are known both for conserved and nonconserved 
cases 0, H| and an intuitive deconfinement of boundary 
layer approach provides a physical picture of the phase 
transitions |f| EtHl ■ Several variants of ASEP have also 
been studied[9l liol. [Tl| . 

For 1-D ASEP chains, the phase diagram for the case 
with conservation in the bulk is shown in Fig[TJb). For 
large length, the phases are characterized by the density 
p{x),x £ [0,1]. The external drives at the boundaries 
maintain a density p = a at x = and p = 7 = 1 — /? 
at x = 1, and determine the fate of the bulk phase. Un- 
like equilibrium situations, the information of the bulk 
phases and phase transitions are contained in the bound- 
ary behaviour. (This can be termed a "holographic prin- 
ciple"). In the a-phase of FigOJb), t ne bulk density is 
p(x) — a with a thin boundary layer maintaining the den- 
sity at the other end. Similarly, in the 7-phase, p(x) = 7 
in the bulk with a boundary layer at the x = end. 
There is a maximal current phase with p{x) = 0.5 for 
a > 0.5,7 ^ 0.5 with boundary layers on each side pro- 
tecting the bulk. In all these cases, the boundary lay- 
ers are attached to the edges. On the first order phase 
boundary between the a- and the 7-phase, the density 
profile is p(x) = a + (7 — a)x without any boundary 
layer. In case of a non-conservation in the bulk, this 
phase boundary gets replaced by a shock phase with 
localized shocks on the trackQ. This additional shock 
phase can be understood as a deconfinement transition 
of the shock from the boundary Q. 

Here, we consider a collection of such one dimensional 
ASEP chains diffusively coupled to form a two dimen- 
sional ASEP (2-D ASEP). The transverse diffusion does 
not lead to any current in the extra dimension but affects 
the bulk and boundary in the preferred forward direc- 
tion. An arbitrary chain may seem to have nonconser- 
vation through the leakage to or from the neighbouring 
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FIG. 2: (a): Spatial density distribution p(i,j) for a particular 
j along the forward direction i for various values of p (a = 
7 = 0.2). The vertical width of any point on the line reflects 
the variation in the density in the transverse direction. For 
p < p c curves form one group and for p > p c the bulk density 
is determined by the left boundary p{l,j) = a/p. (b): The 
presence of diffusing shocks for a — 0.2 and 7 = 0.8 when 

p > Pc- 



chains but there is an overall bulk conservation on the 
lattice. We show here from simulations the existence of 
the maximal current phase with p = 0.5 for high trans- 
verse coupling over a wider range of a and 7 with a phase 
transition to the conventional phase at a critical coupling. 
The phase transition behaviour in this situation can be 
analyzed through the changes in the boundary layers. To 
do so, we also consider a few variants of the model both 
in one and two dimensions. 



II. MODEL 

Consider a modified asymmetric exclusion process 
(ASEP) on a sheet ofuixn sites as shown in Fig. []Ja) 
with forward particle jump probability p > and the 
transverse (perpendicular to the forward direction) prob- 
ability q (with the constraint p + 2q — 1) provided the 
neighboring sites are empty. Here q is a measure of the 
transverse coupling of the chains. For q = 0, we get back 
independent 1-D ASEP chains. On the left boundary 
i = 1, 1 < j ' < n particles are injected at a rate a and 
on the right boundary (i = m, 1 < j < n) particles are 
withdrawn at a rate (3. The sheet is folded in a cylin- 
drical geometry to impose periodic boundary conditions 
in the transverse direction i.e., sites (i,j — 1) are iden- 
tified with sites (i, j = n + 1). Thus in the steady state 
situation we have a net particle current in the forward 
direction only, and no particle current in the transverse 
direction, because the probabilities of up- and down- hops 
are the same. 



A. Mean field analysis 

The occupation number at site (i,j) is t^j = or 1 
depending upon whether the site is empty or occupied. 
The rate equation governing the average particle density 
distribution p(x, y) = (r^j) (where the average is over all 
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FIG. 3: The density phase transition for various values of a 
and (5 in (a) 2-D and (b) 1-D cases. The average bulk density 
p — 0.5 up to p < p c , but varies with p for p > p c . For 
both, curve (1) is for a = 0.8, /3 = 0.2, curve (2) for a = 0.2, 
(3 = 0.2, curve (3) for a = 0.8, (3 = 0.8, and curve (4) is for 
a = 0.2, P = 0.8, as marked in Fig. [j] (b). All these have 
same p c ■ 



realizations of the process) in the bulk is: 
,,T " P [((l njyn-u) ((1 - T i+1J )ru)] + Jl 



(It 



where the transverse part is 

J?i = ?[-((! - n,j+i)n,j) - ((1 - Ti,j-i)T itj ) 



(la) 



+((! - n,j)Ti,j+i) + ((1 - n^nj-! 

The rate equation for the two boundaries are 



>]• (lb) 



dt 



+a(l - nj) -p((l - T 2 , j )r ld ) + J%, (lc) 



= -P(r m ^+P((l-T m , 3 )r m ^) + Jl r (Id) 

It is interesting to note an invariance in the above equa- 
tions known as the particle-hole symmetry. It implies 
that if we change a to 1 — /3 and /3 to 1 — a with r 
changed to 1 — r, the equations of the process remains 
invariant. 

In a mean-field independent-site approximation, one 
sees that p ^constant is a solution of the bulk equation 
in the steady state. The phase of the system is then 
determined by the boundary conditions. It transpires 
that a constant density cannot satisfy in general both 
the boundary conditions. This importance of the bound- 
ary, i.e., the choice of one, both or none of the bound- 
ary conditions, is at the heart of the phase transitions. 
One can in addition do a stability analysis to see that 
a constant bulk density is indeed a stable solution [l3|. 
With the Boltzmann approximation (neglecting near- 
est neighbour correlations) i.e., (^^-(1 — Ti±i t j±i)) = 
fa.iK 1 ~ (n±i,j±i))„ we take (nj) = p + Sp(i,j,t), 
with Sp a small perturbation. In terms of the Fourier 
modes, 



5p(k, t)=J2 e- lk * x ~ zk « y Sp(x, y, t), k x(y) = 



x.y 



(2) 
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where L is the dimension of the lattice, and k denotes 
{k x , k y }, Eq. [la] can be written as 



ddp(k, t) 
Jt 



n( P ,g,k)*p(k,t) J 



(3) 



with 



f2(p, q, k) = ipsin k x (2p — 1) +pco$k x + (1— p) cos k y — 1. 

(4) 

Since pcos fc^ + (1 — p) cos fc y — 1 < 0, the negativity of 
the real part of f2 insures decaying perturbations and 
stability. This linear stability analysis, though useful 
in the context of traffic jams in similar two-dimensional 
models 13], is not enough for ASEP. 



B. Simulation 

To simulate the process for any p we use a random 
sequential update scheme. Starting from a random dis- 
tribution, we allow the system to reach a steady state. 
From the simulations we study the spatial density distri- 
bution and currents for various values of p and for various 
sizes of the lattice. 

To analyze in detail the density dependence on p, let 
us define the average bulk density p(p) (for given a and 
P), 



1 N 1 

cyc—1 i,j(ZA ceri 



' h3 



(p). 



(5) 



The averaging in Eq.(5) is done on a strip (A center ) at 
the center of the cylinder i.e., (m/2 — 4 < i < m/2 + 4, 1 < 
j < n)), and n' (= 9 x n in our case) is the number of sites 
in the central strip A center . N in the above expression is 
the total number of cycles of the simulation (~ 10 6 ) used 
for averaging. 



III. RESULTS 
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FIG. 4: Finite size effects near the transition for the 2-D case 
(a) and for 1-D (b). Data collapse is shown in (c) and (d), 
with f x = (p — p c )L 1/v and f y = (p — p c )L^^. 



For given a and (3, the steady state profiles are of two 
types as shown in Fig. Ela) for several values of p. We see 
that for small p the bulk reaches half-filling and changes 
over to a boundary dependent density for larger p. In 
Fig. [3] p is plotted as a function of p for various values 
of a and j3. The behaviour shown in Fig. [Ha) is evident 
here. For the case a > 0.5, j3 < 0.5, the behaviour is 
complementary to the case a < 0.5, j3 > 0.5. However, 
there is no such transition for a > 0.5, and (3 > 0.5. The 
critical p c also depends upon the values of a and (3 as 
shown in the phase diagram Fig. [ljb). We have studied 
an equivalent 1-D model, because the transverse periodic 
boundary conditions in 2-D has some similarity with 1- 
D. To mimic the behaviour we modify the 1-D ASEP so 
that a particle jumps to the next empty site not with 
probability 1 but with probability p, i.e., particle waits 



with probability 1 — p. The average bulk density in this 
1-D case also shows behaviour similar to the 2-D case and 
is shown in Fig. (3]Jb). One sees the phase transition with 
p. Thus, we have the following three main observations: 

1. p(i,j) — 0.5 for all p less than p c , for all a and (3. 

2. In the regime p > p c , mean-field continuum ap- 
proximation is valid and phase diagram resembles 
the 1-D phase diagram. 

3. In the shaded region marked a-phase in Fig. |Tf b) 

p c = 2a and in the 7-phase region p c = 2/3. 

These results can be explained by examining the bound- 
ary densities. If we do a mean-field approximation in the 
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FIG. 5: (a) Current I(p) vs. p for 2-D case. Curve marked 

(1) is for a = 0.8, /3 = 0.8, curves marked (2) and (3) are 
for a — 0.2, /3 = 0.8 and a — 0.8, /3 — 0.2 respectively, and 
(4) for a = 0.2, /3 = 0.2. Currents are rescaled to match the 
bulk value p(l — p) at p — 1.0, to correct for finite size effects. 
Similar result holds for currents in 1-D case also. The solid 
line is I(p) = pp(l — p) for the parameters of curves (2) and 
(3). In (b) the dependence of p c on a and (5 is shown. For 
all the curves, a + (3 = 1. The curves are: (1) a — 0.85, 

(2) a = 0.15,(3) a = 0.75, (4) a = 0.25, (5) a = 0.65, (6) 
a = 0.35, (7) a = 0.55, (8) a = 0.45, The main observation 
is that p c = 2a for downward curves (a-phase region of Fig. 
QJb)) and p c = 2/3 for the up going curves (7-phase region of 

pig. up))). 



steady-state situation of Eq. [Tel then a homogeneous 
density would give (a — pp){l — p) = or p = a/p for 
the left boundary. The bulk current is expected to be 
I(p) — PP(1 — p) ( see below) as shown in Fig. 02a). 
For ASEP, the bulk satisfies the left boundary condition 
only in the a-phase which requires the boundary density 
to be less than or equal to 0.5. Therefore a maximal cur- 
rent phase is expected if a/p > 0.5, i.e., p c = 2a. The 
left boundary layer then develops (FigG^a)) for p < p c . 
The density variation in the boundary layer vitiates the 
simple argument because the density gradient dependent 
diffusive part of the boundary current needs to be taken 
into account. The net boundary density is obtained by 
the balance of the input and the outflow consisting of the 
hopping and the diffusive parts. Similar argument holds 
in the 7-phase region for and p c . For the 7-phase, the 
right density is p = [7— (1 — p)]/p if there is no boundary 
layer. The bulk density is controlled by this boundary 
value (rather than the withdrawal rate) so that it also 
takes the same value as the boundary. These observa- 
tions are supported by Fig. E2a,b). 

It is known for ASEP, that on the first order phase 
boundary separating the a- and the 7-phases, there are 
shocks that diffuse slowly on the track vanishing or get- 
ting created at the boundaries only. Same thing happens 
here also on the phase boundary which is still set by 
a = l — 7. Because of slow diffusion of the shock, the 
measured density in the central patch could be either 
that of the a-phase or of the 7-phase. This is shown in 
Figs. (a) and (b)). The density remains constant for 
p < p c , but after this (p > p c ) average density shows 
an erratic behaviour, fluctuating wildly [|. The special 
point where the three phase boundaries meet is now at 
a = p/2, 7=1 — (p/2) in the a, 7 plane. 



The above meanfield results seem to suggest a singu- 
larity in the density as a function of p, because p = 0.5 
for p < p c but p = a/p for p > p c . Such a singularity 
is expected only in the long chain limit (infinitely long 
system) and not in finite systems. Fig. |4ja,b) shows a 
strong size dependence near p c . For equilibrium phase 
transitions, singularities are rounded off by finite size 
when the size of the system is comparable to the char- 
acteristic length scale for the transition. The finite size 
behaviour, especially the size dependence, then follows a 
finite size scaling form. In that spirit, let us make a finite 
size scaling ansatz for this nonequilibrium case as 

p-pc-L-vfdp-p^L 1 /"), (6) 

with, 

p-pc~\p-p c r forL^oc (7) 

where p c = 0.5 is the constant density for p < p c , L is the 
linear dimension of the system (lattice or chain), and, p, 
and v are scaling indices, then, to recover the meanfield 
results, we need to have p = v~ x . We have used the 
Bhattacharjee-Seno method for data-collapse (l2j . In Fig. 
UJc) the data collapse scaling is shown for 2-D for which 
we get pc = 0.40 ± 0.008, p = 0.69 ± 0.07, V 1 = 
0.72 ± 0.03. For the 1-D case (Fig. gfd)), we have p c = 
0.401 ± 0.006, p = 0.46 ± 0.02, u' 1 = 0.44 ± 0.06. 
These are consistent with the prediction of pv = 1 . The 
characteristic length scale seems to diverge as £ ~| p — 
Pc \~ v which is set by the width of the boundary layer. 
Meanfield analysis is not fine enough to get this length 
properly. 

Since the current is a measure of jumps from occu- 
pied sites to nearest vacant site in the forward direction, 
the probability of site occupation is p, the probability 
of vacancy of the next site is 1 — p, and jump probabil- 
ity in the forward direction is p, thus, the net current in 
the forward direction is I(p) = pp(l — p). Consequently, 
I(p) = p/A for p < pc, while I(p) = a(l — a/p) for p > p c , 
joining continuously at p = p c with a slope discontinuity. 
Fig. GHa) shows the overall agreement of the measured 
current and this general form of the current when the 
correnponding p obtained from the simulation is used. 
However, finite size rounding masks the expected singu- 
larity at p — p c in this current plot. 

In order to show that the above results, though bound- 
ary driven, are not a consequence of local perturbations 
at the boundary, we considered a variant of the model 
where the transverse coupling is only at the two ends. 
We have put p — 1 in the all the bulk sites i.e., for sites 
2 < i < n — 1 and kept finite p for the first and the 
last site, i.e., jumps from first site to second and n — 1th 
site to nth happen with finite p. We see that the sys- 
tem self-organizes to a state with new boundaries that 
control the bulk density. The actual drives (the injec- 
tion and withdrawal rates) passively help in creating the 
relevant boundary conditions. In particular, we observe 
that the transition induced by p for the bulk case is no 
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FIG. 6: The importance of bulk and boundary (a) the density 
transition is missing when we consider only boundary sites 
(finite p) and make p = 1 in all interior sites, (b) spatial 
density profiles in 2-D case for various values of p. For a, j = 
0.2. Similar behaviour is observed in equivalent 1-D case. 



longer present. The behavior of average p with p is shown 
in Fig. [6] (a) and the corresponding density profiles are 
shown in Fig. [f^b) (similar profile has benn observed in 1- 
D case also) . The behaviour of average p with very small 
p < 0.05 shows a long living transient state, due to the 
very small forward motion. These observations indicate 
that the transition is due to a co-operative phenomenon, 
where bulk and boundary play their role co-operatively 
and inter-dependent way. 



IV. SUMMARY 



In conclusion, the continuous transition from the in- 
jection rate dominated phase to the maximal current 
phase has been observed as a function of forward transi- 
tion probability p in a two dimensional ASEP (diffusively 
coupled chains). The transition shows finite size effects, 
reminiscent of equilibrium phase transitions, and finite 
size scaling predicts exponents which are consistent with 
the mean field theory predictions. The bottleneck cre- 
ated at the boundary by the transverse coupling changes 
the effective particle densities at the two boundaries and 
the ensuing phase diagram can then be mapped out from 
the 1-d phase diagrams with p = 1, with the multicrit- 
ical point shifting to (cv = |,7 = 1 — |). However no 
such transition can be induced if artificial bottlenecks 
are created at the boundaries only. In such situations, 
the particles organize themselves to form a new or effec- 
tive boundary density which then as per the holographic 
principle fixes the bulk density. This reiterates that the 
nonequilibrium transitions observed are cooperative but 
boundary driven and the boundary layers contain the in- 
formation about the bulk. 
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